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We present a complete, next-to- leading-order (NLO), 
leading-twist QCD analysis of deeply virtual Compton scat- 
tering (DVCS) observables, in the MS scheme, and in the 
kinematic ranges of the HI, ZEUS and HERMES experi- 
ments. We use a modified form of Radyushkin's ansatz for 
the input model for the generalized parton distributions. We 
present results for leading order (LO) and NLO for represen- 
tative observables and find that they compare favourably to 
the available data. 

PACS numbers: ll.lO.Hi, 12.38.Bx, 13.60.-r 

Deeply virtual Compton scattering (DVCS) 
7*('Z) + p{P) 7(9') + p{P')i is the most promising 
process for accessing generalized parton distributions 
(GPDs) |l]J^,|^,|lO|"p^ which carry new information about 
the dynamical degrees of freedom inside a nucleon. GPDs 
are an extension of the well-known parton distribution 
functions (PDFs) appearing in inclusive processes and 
are defined as the Fourier transform of nonlocal light- 
cone operators sandwiched between nucleon states of dif- 
ferent momenta, commensurate with a finite momen- 
tum transfer in the t-channel. These distributions are 
true two-particle correlation functions and contain, in 
addition to the usual PDF-type information residing in 
the so-called "Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) region" supplementary information about 
the distribution amplitudes of virtual "meson-like" states 
in the nucleon in the so-called "Efremov-Radyushkin- 
Brodsky-Lepage (ERBL) region" We have recently 
presented a full numerical solution of the associated 
renormalization group equations at next-to-leading order 
(NLO) accuracy, for both the polarized and unpolarized 
distributions |15| . 

In this paper, we report on a full NLO, leading- twist 
QCD analysis carried out for unpolarized and singly po- 
larized (with respect to the lepton) DVCS and its phys- 
ical observables based on the DVCS factorization the- 
orem [0,1). The precise technical details of this analy- 
sis are presented in considerably more detail elsewhere 
psj-pTt- DVCS and its observables are accessed experi- 
mentally in deep inelastic scattering experiments B us- 



ing the leptonic process e p 



e^pj. 



Hence, on the lep- 



ton level, DVCS interferes with the Bethe-Heitler (BH) 
process which has the same final state (for BH, the fi- 



nal state photon is radiated off either the initial or final 
state lepton). The exact DVCS kinematics we base our 
analysis on can be found in [p^ . The physical observ- 
ables we use as examples are the single spin asymmetry 
(SSA) and the (unpolarized) azimuthal angle asymmetry 
(AAA) as well as the photon level cross section, which 
we will define below. We will demonstrate the feasibility 
of constraining the GPDs from experimental data using 
these observables, as well as the relative size of the NLO 
corrections to the observables. 

For our input GPDs we use a modified form of the dou- 
ble distribution ansatz due to Radyushkin which 
assumes that the dependence on t = {P ~ P'Y fac- 
torizes. In ]l7[ | we demonstrated that the unmodified 
ansatz overshoots the DESY ep collider HERA data, es- 
sentially because it samples an extrapolation of the in- 
clusive unpolarized PDFs, which rise steeply as x de- 
creases, right down to a; = 0. For the quark distribution 
q{x) = xq{x)/x ^ x^^''^ , this leads to a large enhance- 
ment of the quark GPD relative to the PDF at the input 
scale. Our modification is to introduce a lower cutoff 
on X, which we motivate below, and it leads to a good 
description of the available data. Like inclusive PDFs, 
the input GPDs are unknown nonperturbative bound- 
ary conditions which should ultimately be constrained 
by a QCD analysis of the data. The HI data already 
appears to be ruling out Radyushkin's ansatz in its un- 
modified form. NLO analyses of the evolution of GPDs 
and the DVCS amplitudes themselves have been carried 
out previously for large Xbj 0,|o| and at LO for DVCS 
observables ||7|,|2l|. 

The GPD combinations, appearing in the representa- 
tion of the factorization theorem j^,^ for the DVCS am- 
plitude used in our analysis, are the singlet-type combi- 
nations for each parton species, which are defined by 

:FSia),V/A^X,C,ll',t)^ 



(l-C/2) 
.F«^^/^(X,C,M',i) = 



(1) 
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(l-C/2) 

where ( is the skewedness parameter, v — {X ~C/2) /(I — 
C/2), fj, the factorization scale, a = u,d, s labels the quark 
flavor, and V stands for the unpolarized (vector) case and 
A for the polarized (axial-vector) case. 

Our representation of the GPDs is identical to the non- 
diagonal representation defined in p^] and is different 
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from the usual one (see for example ||Tl[), which is de- 
fined symmetrically with respect to the incoming and 
outgoing nucleon plus momentum (defined on the inter- 
val V S J— 1, 1] and symmetric about u = 0). The GPDs 
in Eq. (Q) have plus momentum fractions (on the inter- 
val X G [0, 1]) with respect to the incoming nucleon mo- 
mentum, P, in analogy to the PDFs of inclusive reac- 
tions, with the ERBL region in the interval X £ [0, C] 
and the DGLAP region in the interval X G [C, 1]. The 
transformation between the symmetric and nondiagonal 
representation is given in p^ ]. Furthermore, within the 
nondiagonal representation C = Xi,j = —q^/2P ■ q. The 
symmetries of the GPDs, which were previously man- 
ifest about V = 0, are now manifest about the point 
X = C/2. Note that J^'^-^ and T^^"^ are therefore an- 
tisymmetric and J-^'^ and J^^'^ are symmetric about 
the point X = C/2. 

The DVCS amplitude, Tuvcs{l*P ~* ip)i factorizes 
into a convolution of GPDs with coefhcient func- 
tions: 



c 



P.V. f dX T^^^/-^ (2X/C - 1) Jc-5(")'^/^(X, C, Q^ t) 
Jo 



T j dXT^'^/^(l-2X/C)^^(")'''/^(X,C,Q',t) 
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DVCS 
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2-C 

c 



P.V. f dX T<^'^/^ {2X/C - I) T3'^/'^{XX,Q^,t) 
Jo 



± J dX PS'^/^ {l-2X/C)TS'^^^{XX,Q^,t) 



(2) 



with the choice fj? = Q'^ = —q^ and the T's are the LO or 
NLO coefficient functions taken from Eqs. (14) - (17) of 
pot - P.V. stands for the Cauchy principal value prescrip- 
tion, which implements the -fie prescription to regulate 
the divergences in T{2X/C, — 1) at X = ^. Depending 
on the region of integration the coefficient functions have 
real and imaginary parts, which in turn generate real and 
imaginary parts of the DVCS amplitudes. 

The triple differential cross section, on the lepton level, 
contains pure BH and DVCS terms, and an interference 
term. It is given in general by 



da(3)(e±p 



ot' 



dxhjdQ'^ d\t\ 



27r 



dxhidQ"^ d\t\d<^ 



-1/2 



where 

ir±P : 



DVCS\ 



±(r, 



DVCs'^BH 



T, 



DVCS'^BH 



±|2 



(3) 



and (j> is the relative angle between the lepton and hadron 
scattering planes (we use the special Lorentz frame of [T^ ] 
to define 4>), y — Q^/x^jS is the energy fraction lost by 
the incoming lepton, S is the total center of mass energy 
and M is the nucleon mass. The exact, leading twist, 
expressions for the DVCS square, interference and BH 
square, both for an unpolarized and longitudinally po- 
larized probe on an unpolarized target, were taken from 
Eqs. (24) - (32) of lljl with the full dependence of the 
BH propagators restored [ pTj . 

The (unpolarized) azimuthal angle asymmetry (AAA) 
and the single spin asymmetry (SSA) are defined by 



AAA 



and 



SSA 



Jo^''d(/)sin(/)AaW 



(4) 



(5) 



where Acr^''-' = da^ ~da^ with | and | signifying that the 
electron or positron beam is polarized along or against 
the beam direction, respectively. These definitions make 
the asymmetries directly proportional to the real part 
of a combination of DVCS amplitudes, for AAA, and the 
imaginary part of a combination of DVCS amplitudes, for 
SSA. For small x and t, these combinations of amplitudes 
reduce to just the unpolarized helicity nonflip amplitude 
(which are obtained from a form factor decomposition of 
the H functions defined above) for both the AAA and 
SSA [|l]j2|. 

As a model for our input GPDs we use the double dis- 
tribution (DD) representation |l| with a factorized ansatz 
for the t-dependence |l^] of the DDs, which has been the 
model chosen for virtually all phenomenological studies 
of DVCS: 



Fl(^y^'''/^{x',y',Qlt) = 

TT^^'ix', y') p'^<')la-^l^{x\ Ql) r'^/ay/^it) 



(6) 



BH\ 



where the profile functions tt'^'^{x' ,y') are asymptotic 
shape functions |Q| for quarks and gluons. The 
jq{a)/g,v/A g^j.g unpolarized/polarized inclusive input 
PDFs defined at the input scale Qq. The factor- 
ized t-dependencc of the DDs is given by form factors, 
ri-3y/^{t)^ depending on the parton species and whether 
the proton flips its helicity or not p]. Here, we follow 
the model specified in Section. (4) of |21|], with a further 
assumption that the gluon distributions have the same 
^-dependence as the sea |l^. The GPDs that we use for 
evolution are the nondiagonal Ts of Eq.(|^): thus to de- 
fine the input model we require the reduction formula 
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relating the ff(w)-functions to the DDs: 



HivX) 



1 /-i-b' 
dx' / dy' 

-1 J-l + \x'\ 



S{x' + 



2-C 



V Fdd{x ,y). (7) 



Performing the y' integration, using the d function, mod- 
ifies the hmits on the x integral according to the region. 
For example, in the DGLAP region, X > (, one has for 
the quark distribution 



l-C 



c 



(8) 



so that a,s X ^ ( the region x' ^ is probed in the 
PDF. For the polarized case this presents no problem 
since the available polarized distributions tend to zero 
as x' — > and this region is irrelevant to the integral. 
However, small x inclusive unpolarized PDFs are known 
to increase rapidly as x decreases so that this model is 
very sensitive to the extrapolation of the PDFs into the 
unmeasured region. Taking the model at face value leads 
to quark GPDs which are much larger than PDFs close 
to X = C and overshoots the HERA data on DVCS 
This region is also probed in the ERBL region. 

However, if one considers the imaginary parts of par- 
ticular DVCS graphs (which correspond to quarks and 
gluons on mass shell in the final state), one realizes 
that the limit X cannot be realized in practice 

because the partons are obliged to produce finite mass 
hadrons in the final state. For example, for the pro- 
qq, we have {XP + q)^ 



cess f/7 



M^^ which at 



the input scale, Qo, implies X — ( > M^^C/Qo- 
fact one can argue that X ~ C, > for each of the 
graphs, where a ~ ^^("t^orfron/Qo)- other words, fi- 
nite mass effects become important for X = C. As a 
phenomenological model, we decide to pick the generic 
hadronic mass to equal the lightest vector meson mass, 
nT-hadron = "mp — 0.77 GcV, which gives a ~ 1/2 for 
Qo = 1 GeV, and replace the lower limit of the inte- 
gral in Eq.(|^) (and all similar integrals) with aC ~ C/2- 
This change [g4| ensures that the PDFs are only sam- 
pled at x' ^ C(C) and above and thereby reduces the 
enhancement of the GPDs relative to the PDFs (see ||2^ ] 
for further discussion). 

In addition to these contributions from the double dis- 
tributions the unpolarized singlet GPDs also contain a 
"D term" ^j^, which is only nonzero in the ERBL 
region, and ensures the correct polynomiality in C, p^ ] 
of the GPD and has a physical interpretation as a vir- 
tual meson exchange in the t channel. The unpolarized 
helicity-fiip GPD, due to lack of knowledge, is assumed 
to contain only the D term with an overall minus sign 
as required by the first moment sum rule for the sum of 



unpolarized helicity-flip and nonflip, in which the indi- 
vidual D terms must cancel. The polarized helicity-fiip 
GPD was taken to be the asymptotic pion distribution 
amplitude due to the presence of the pion pole in this 
GPD, thus it does not evolve. 

Within this class of input model, we specify a particu- 
lar input model for the GPDs by using the latest Martin- 
Roberts-Stirling-Thorne (MRST) LO and NLO sets g] 
for the unpolarized PDFs (which have been tuned to the 
latest inclusive data from HERA) and the Gliick-Reya- 
Stratman- Vogelsang (GRSV) | 29[ set for the polarized 
partons (for use in Eqs. (|l|) (|) (|7|) ) at the input scale 
Qq = 1 GeV^. We illustrated the spread of predictions 
one finds using different input PDFs in [0. In the cal- 
culation of these observables, only gluons, up, down and 
strange quarks were used. The small charm quark con- 
tribution, generated by evolution, has justifiably been 
neglected so far (see [|l^ for more details). 

The input GPDs are evolved in both LO and NLO us- 
ing an evolution program based on direct numerical in- 
tegration, specifically written for this purpose (l^. The 
DVCS amplitudes are then computed via numerical inte- 
gration, using LO or NLO coefficient functions in Eq.(||) 
as appropriate and we chose the renormalization scale to 
be equal to the factorization scale, = /i^- It turns out 
that the theoretical uncertainties due to scale variation 
are smaller than the uncertainties due to different inputs, 
therefore we will give our results without error bands for 
the scale dependence (see for more details). The 

resultant DVCS amplitudes, calculated to both LO and 
NLO accuracy, using the respective LO and NLO input 
models, were the input for a further numerical program 
which computes the final DVCS observables (for details 
see ) . We will publish all of these codes on the Inter- 
net for use by the community jsotl . 

In Fig. 1^ we show representative NLO results for 
the azimuthal angle and single spin asymmetries, which 
by construction are directly sensitive to linear combi- 
nations of real and imaginary parts of DVCS ampli- 
tudes, for two values of the skewedness parameter, i.e. 
C = Xbj — 0.2,0.0002, which are indicative of values ac- 
cessible at HERMES and ZEUS/Hl, respectively. We 
note that the NLO predictions for both observables are 
discernibly different from the LO ones, although the NLO 
corrections to the single spin asymmetry are never more 
than about 20%, whereas the NLO corrections to the 
azimuthal angle asymmetry can be larger than 100%. 
This shows that the appearance of the gluon at NLO 
is stronger in the real part than in the imaginary part of 
the amplitude (see also ]l6t). 

The HERMES Collaboration ||] quotes two values for 
SSA: SSA = -0.18 ± 0.05 ± 0.05 and (SSA) = -0.23 ± 
0.03 ± 0.04. The first, which assumes that their missing 
mass equals the proton mass, is quoted at the following 
average values: < x >= 0.11, < >= 2.5 GeV^ and 
< t >= —0.27 GeV^. The second (average) value is the 
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SSA integrated over the missing mass at the same average 
values of and t. At this point in x,Q^,t, we find 

SSA = -0.31 in LO and SSA = -0.20 in NLO for 
MRSTOl/GRSVOO input PDFs, i.e. good agreement. 
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FIG. 1. The single spin and azimuthal angle asymmetries 
as functions of t for ZEUS/Hl (C = Xbj = 0.0002) and HER- 
MES (C = Xbj = 0.2 ) kinematics, using MRSTOl (unpolar- 
ized) and GRSV (polarized) input PDFs. The solid curves 
are the full LO results and the dashed curves show the full 
NLO results for the respective input GPDs. 

We now compare our model to the recent HI data on 
the cj{'j*p — > 7p) photon level cross section in Fig. ||, 
defined through 



LO 




NLO 




H1 data -+ : 















60 80 1 00 

W (GeV) 




10 12 

Q"2 (GeV"2) 

FIG. 2. Photon level cross section a{"/''P), at 
Q2 = 4.5 GeV^ in W (upper figure) and at VF = 75 GeV 
in Q2 (lower figure), for MRST/GRSV at LO and NLO. 

We have modified Radyushkin's ansatz for the GPDs 
by introducing a ^-dependent lower cutoff, « (^/2, on 
the sampling of the PDFs at small x. Such a cutoff may 
be justified by arguing that the region X ^ C is sup- 
pressed by finite mass effects. Without further tuning of 
a ^ 0{mf^^j^j.g^/ Qq) , we can produce a good description 
of the current data using a NLO QCD model. As the 
precision of the data improves it may become necessary 
to use this cutoff as a fit parameter to constrain the in- 
put GPDs. We compare our NLO model to the LO result 
to quantify the size of the NLO QCD corrections to the 
DVCS observables. 

We thank M. Diehl, D. Midler and C. Weiss for help- 
ful discussions. A.F. was supported by the DFG (FR 
1524/1-1), and M.M. by PPARC. 



\T, 



DVCS\ t=0, 



(9) 



where B stems from the t integration and, within our 
model for the t dependence, was found to be about 
B w 6.5 GeV^^. With our modification to Radyushkin's 
ansatz both LO and NLO curves appear to agree well 
with the data. Note the decrease of about 20 — 40% in 
going from LO to NLO, which is due to a large, negative 
gluon amplitude and in line with results from inclusive 
reactions. 
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